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ABSTRACT 

The  experimental  data  on  temperature,  composition, 
and  aerodynamic  profiles  through  a  flat,  premixed  CH^-O^ 
flame  obtained  by  techniques  described  in  the  report  CM-968, 
Characteristic  Profiles  in  a  Low-Pressure ,  Laminar ,  Lean, 
Premixed  Me thane -Oxygen  Flame ,  of  this  series  are  analyzed 
in  terms  of  matter  and  energy  conservation.  The  composi¬ 
tion  profiles  for  the  seven  stable  species  found  in  the 
flame  §.re  corrected  for  molecular  diffusion,  using  dif¬ 
fusion  coefficients  measured  in  this  Laboratory,  and  the 
resulting  flux  distributions  are  shown  to  give  a  satis¬ 
factory  balance  of  carbon  and  hydrogen  through  the  flame. 

The  fluxes  of  enthalpy  by  means  of  convection,  diffusion, 
and  conduction  are  given  separately  and  summed  to  show 
the  degree  to  which  energy  conservation  is  fulfilled  by 
the  data.  The  various  approximations  used  in  the  treat¬ 
ment  are  critically  discussed,  and  it  is  shown  that  the 
general  reliability  of  the  results  has  been  considerably 
improved  over  previous  laminar  flame  structure  analyses. 
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METHANE-OXYGEN  FLAME  STRUCTURE 
CONSERVATION  OF  MATTER  AND  ENERGY  IN  THE  ONE-TENTH 
ATMOSPHERE  FLAME 

I .  INTRODUCTION 

The  development  of  experimental  techniques  for  the 
measurement  of  temperature,  composition,  and  aerodynamic 
profiles  through  laminar  flame  zones  has  been  reported  in 
earlier  publications  of  a  series  from  the  Applied  Physics 
Laboratory  (Refs.  1,  2  and  3).  This  work  culminated  in  a 
set  of  data  obtained  in  a  stoichiometric,  premixed,  pro¬ 
pane-air  flame  of  conical  geometry  at  1/4  atmosphere  pres¬ 
sure.  The  fourth  paper  of  the  Combustion  and  Flame  series 
(Ref.  4)  reported  the  analysis  of  these  data  in  terms  of 
material  transport. 

The  experience  gained  up  to  that  time  pointed  the 
way  to  various  refinements  and  changes  in  the  experimental 
techniques,  and  also  to  the  desirability  of  studying  a 
flame  system  which  was  chemically  simpler  insofar  as  the 
number  of  stable  species  present  is  concerned.  The  changes 
in  technique  were  described  in  detail  in  Ref.  5.  In  brief 
summary,  these  included  the  following: 

1.  A  flat-flame  stabilized  on  a  screen  burner  was 
used  instead  of  the  conical  flame  used  in  the 
previous  work.  This  brought  about  some  simpli¬ 
fication  in  the  geometrical  handling  of  the  data, 
and  in  the  physical  manipulation  of  the  apparatus. 
References  may  be  found  on  pages  23,  24, 'and  25. 
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2.  The  operating  pressure  was  lower  (0.1  atmosphere 
instead  of  0.25  atmosphere)  which  gave  a  thicker 
flame  and  hence  eased  the  requirements  on  spatial 
resolution . 

3.  Gas  sampling  and  mass  spectrometer  analyses  were 
done  in  a  continuous  flow  system  instead  of  batch- 
wise.  This  afforded  much  better  precision,  and 
circumvented  wall  adsorption  difficulties.  In 
particular,  this  system  permitted  good  analyses 
for  water  to  be  realized. 

4.  Temperature  profiles  were  obtained  with  a  fine, 
silica-coated  Pt-Pt,  10  per  cent  Rh  thermocouple 
rather  than  by  particle-track  photography  (which 
was  used  only  as  a  supplementary  technique).  Con¬ 
siderably  better  precision  was  thus  attained. 

5.  Improvements  in  controlling  and  monitoring  the 
flow  of  gases  to  the  burner  were  introduced  which 
led  to  better  flame  stability.  A  monitoring  ther¬ 
mocouple  was  mounted  just  downstream  of  the  screen, 
which  allowed  the  flame  position  to  be  monitored 

at  all  times. 

The  chemically  simpler  flame  system  finally  decided 
on  was  a  mixture  of  methane  and  oxygen.  It  was  thought  that 
this  represented  a  reasonable  compromise  between  chemical 
simplicity,  ease  of  handling,  and  some  degree  of  general 
interest  as  a  common  combustible  mixture.  As  shown  in  Ref . 

5,  this  flame  contained  five  major  stable  species  (methane, 
oxygen,  carbon  monoxide,  carbon  dioxide,  and  water)  and  two 
minor  ones  (hydrogen  and  formaldehyde)  for  which  analyses 
were  made,  which  is  considerably  simpler  than  the  twelve 
components  found  in  the  propane-air  flame.  No  other  species 
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(other  than  impurities)  at  a  concentration  greater  than 
0.01  per  cent  were  found. 

One  other  important  improvement  incorporated  in  this 
new  work  should  be  mentioned.  This  is  that  experimentally 
determined  molecular  diffusion  coefficients  were  used  in  the 
data  analysis  to  be  described.  In  view  of  the  pronounced 
(in  some  cases  the  predominant)  effect  of  diffusion  on  the 
flux  profiles  of  the  various  flame  constituents,  this  is  a 
decided  gain  in  the  reliability  of  this  general  approach  to 
laminar  flame  studies,  particularly  in  the  derivation  of 
chemical  kinetic  information  from  the  data.  In  the  previous 
work  on  propane-air  flames,  it  was  necessary  to  use  values 
of  the  molecular  diffusion  coefficients  calculated  from 
kinetic  theory  and  empirical  low  temperature  viscosity 
parameters,  so  that  the  reliability  of  £he  diffusion  coef¬ 
ficients  at  flame  temperatures  was  quite  uncertain.  Since 
then  a  new  method  of  measuring  diffusion  coefficients  up 
to  moderately  high  temperatures  (~'1200°K)  has  been  devel¬ 
oped  at  the  ..i;'lied  Physics  Laboratory  (Refs.  6,  7  and  8) 
and  used  specifically  to  measure  the  pertinent  coeffi¬ 
cients  for  the  methane-oxygen  flame  (Ref.  9). 

The  present  paper  describes  the  analysis  of  the  data 
on  the  methane-oxygen  flame  from  the  viewpoint  of  the  trans¬ 
fer  of  mass  and  energy.  It  is  believed  that  the  various 
changes  and  refinements  noted  above  have  contributed  to 
making  the  analysis  of  this  flame  considerably  more  reliable 
than  was  possible  before. 
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II.  PREPARATION  OF  THE  DATA 

As  described  in  detail  in  Ref.  5,  the  first  paper 
of  this  series,  the  flame  used  in  this  study  was  flat, 
and  held  at  a  pressure  of  7.60  cm  Hg.  The  initial  mix¬ 
ture  composition  (mole  fractions)  aside  from  traces  of 
impurities  totaling  less  than  1  per  cent  (argon,  nitro¬ 
gen,  carbon  dioxide,  and  water)  was  as  follows: 

”  0.0785  ;  xA  =  0.9143. 

CH4  °2J 

...  v‘/l 

The  data  finally  used  in  the  subsequent  analysis  con¬ 
sisted  of  a  set  of  concentration  profiles  for  the  stable 
species  present,  a  thermocouple  temperature  profile,  a 
streamtube  area  ratio  profile,  and  a  mass  average  gas 
velocity  profile.  The  latter  was  obtained  from  a  meas¬ 
urement  of  the  final  (hot  boundary)  gas  velocity  vf  by 
means  of  particle-track  photography  and  the  temperature 
profile,  making  use  of  the  continuity  relation 

Pf  vf  af  =  pva  ,  (1) 

where  p  is  density  and  a  is  the  cross-sectional  area  of  a 
streamtube,  and  the  equation  of  state 

p  =  PM/RT  ,  (2) 

where  M  is  the  mean  molecular  weight  and  the  other  symbols 
are  conventional.  Since  the  thermocouple  temperature  pro¬ 
file  was  much  smoother  and  more  reproducible  than  the 
particle-track  velocity  profile,  it  was  felt  that  this  pro¬ 
cedure  of  deriving  the  velocity  from  the  temperature  by 
fitting  at  the  measured  point  where  the  particle-track  data 
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were  most  reliable  was  preferable.  Note  that  this  is  the 
reverse  procedure  to  that  used  previously  (Ref.  4)  where 
the  temperature  profile  was  derived  from  particle-track 
velocity  measurements. 

The  problem  then  arises  as  to  how  best  to  super¬ 
impose,  spatially,  the  temperature  profile  (as  discussed 
above,  the  gas  velocity  and  area  ratio  profiles  were 
coupled  to  the  temperature,  and  thus  not  independently 
variable)  and  the  composition  profiles  determined  by  in¬ 
dependent  means.  It  has  been  pointed  out  (Ref.  1)  that 
when  temperature  and  composition  are  determined  with 
different  devices  (thermocouple  and  sampling  probe)  it 
leads  to  an  uncertainty  in  referring  both  to  a  common 
spatial  origin  even  though  the  positions  of  both  are 
measured  relative  to  a  fixed  point  on  the  burner  (the  moni¬ 
toring  thermocouple) .  The  pneumatic  (sonic  flow)  probe 
temperature  measurements  were  taken  with  the  same  probe  as 
used  in  sampling,  consequently  it  was  assumed  that  a  tem¬ 
perature  read  in  this  way  would  coincide  automatically  with 
the  appropriate  composition  point.  Thus  the  method  used  to 
determine  the  absolute  origin  was  to  superimpose  the  ther¬ 
mocouple  and  pneumatic  probe  temperature  "knees",  that  is, 
the  point  of  intersection  of  straight  lines  extrapolated 
from  the  hot  gas  region  and  the  region  of  rapid  rise  in  the 
temperature  profiles.  The  composition  profile  was  then  fixed 
relative  to  the  thermocouple  temperature  profile.  Additional 
discussion  of  this  problem  is  included  in  the  next  section. 
The  distance  coordinate  z  was  taken  to  be  zero  at  the  moni¬ 
toring  thermocouple  just  downstream  of  the  screen. 
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III.  ANALYSIS  QF  THE  DATA 

Conservation  of  Matter 


With  the  primary  data  on  composition,  temperature, 
and  mass  average  gas  velocity  available,  the  first  step 
in  the  analysis  was  to  take  account  of  diffusion  effects. 
The  quantity  is  the  fraction  of  the  total  mass  flow 
at  any  point  in  the  flame  which  is  due  to  chemical  species 
i,  and  is  given  by 

G1  =  f.Cv+V^/v  ,  (3) 

where  ^  is  the  concentration  expressed  as  mass  fraction 
at  any  point  and  is  the  diffusion  velocity.  It  is  a 
net  flux  variable  which  includes  a  positive  contribution 
from  diffusion  if  the';  concentration  of  i  decreases  down¬ 
stream,  or  a  negative  contribution  in  the  opposite  case. 

As  in  the  previous  work  (Ref.  4),  the  diffusion  velocity 
V\  was  calculated  by  assuming  first,  that  each  species  i 
could  be  regarded  as  a  titace  component  in  the  mixture  and, 
Second,  (since  oxygen  is  present  in  large  excess  everywhere) 
each  species  i  could  be  treated  as  being  in  a  binary  mix¬ 
ture  with  oxygen  (see  Appendix  A).  Under  these  conditions 
the  approximate  relation 

V\  =  -  (Dj/x^  (dx^dz)  (4) 

may  be  used,  where  D.  is  the  binary  diffusion  coefficient 
with  oxygen.  (This  also  neglects  the  effect  of  thermal 
diffusion,  a  point  which  is  discussed  in  Appendix  B.) 
Equation  (4)  was  used  to  compute  the  for  all  the  flame 
species  (except  oxygen,  of  course)  at  each  point  throughout 
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the  flame.  The  diffusion  velocity  for  oxygen  was  then 
obtained  from  the  required  normalization  condition 

2  N.  M.  V.  =  0  (5) 

~  l  l  l 

is  the  molar  density  of  species  i,  and  its  molec¬ 
ular  weight)  which  must  be  fulfilled  in  view  of  the  defi¬ 
nition  of  diffusion  velocity  (relative  to  the  mass  aver¬ 
age  gas  velocity  v). 

The  concentration  gradients  required  in  Eq.  (4)  for 
each  species. were  obtained  by  numerical  differentiation  of 
the  smooth  curves  drawn  through  the  experimental  data,  such 
as  are  shown  as  the  solid  curves  in  Figs.  1  through  7.  The 
numerical  process  was  supplemented  by  graphical  differen¬ 
tiation  at  those  portions  of  the  curve  where  the  slope 
changed  rapidly.  The  internal  consistency  of  the  differen¬ 
tiations  .was  checked.  using  the  requirement  that 
5!  •dx./dz)  =  0,  since  £  x.  =1. 

i  1  i  1  V' 

The  diffusion  coefficients  Di  used  in  Eq.(4)  were 

measured  as  a  function  of  temperature  in  this  Laboratory 
(Ref. 9).  The  point-source  technique  used  permitted  meas¬ 
urement  of  the  gas  pairs  COg-Og  and  HgO-O^  up  to  about 
1100°K,  while  C0~02.  H2~02,  and  CH4-02  were  limited  to 
slightly  lower  temperatures  because  of  self-ignition.  In 
all  cases,  however,  the  temperature  range  covered  was 
enpugh  to  permit  fitting  the  data  to  the  appropriate  kinetic 
theory  formulas  (Ref.  9),  and  reliable  extrapolation  to  the 
highest  temperatures  needed  in  this  flame  (2000°K) .  Dif¬ 
fusion  coefficients  for  the  CH20-02  pair  were  not  measured. 
For  this  minor  component,  the  necessary  values  were  ob¬ 
tained  from  kinetic  theory  (Ref.  10)  using  the  Lennard- 
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Fig.  2  CONCENTRATION  (f)  AND  FLUX  FRACTION  (G) 
PROFILES  OF  WATER  THROUGH  FLAME  ZONE 
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Fig.  3  CONCENTRATION  (f)  AND  FLUX  FRACTION  (G) 
PROFILES  OF  CARBON  DIOXIDE  THROUGH  FLAME  ZONE 
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Fig.  4  CONCENTRATION  (f)  AND  FLUX  FRACTION  (G) 
PROFILES  OF  CARBON  MONOXIDE  THROUGH  FLAME  ZONE 
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Fig.  5  CONCENTRATION  (f)  AND  FLUX  FRACTION  (G) 
PROFILES  OF  HYDROGEN  THROUGH  FLAME  ZONE 
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Fig.  6  CONCENTRATION  (f)  AND  FLUX  FRACTION  (G) 
PROFILES  OF  FORMALDEHYDE  THROUGH  FLAME  ZONE 


TEMPERATURE  (°K) 


Fig.  7  CONCENTRATION  (f)  AND  FLUX  FRACTION  (G) 
PROFILES  OF  OXYGEN  THROUGH  FLAME  ZONE 
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Jones  (6-12)  potential  and  viscosity  parameters.  Since  the 
latter  are  unavailable  for  formaldehyde,  the  values  for 
methanol  were  used,  this  being  a  closely  related  molecule 
of  similar  size  and  mass.  Table  I  lists  the  values  of 
versus  temperature.  The  underlined  values  were  extra¬ 
polated  using  the  Lennard-Jones  (6-12)  potential  with 
parameters  fitted  to  the  lower  temperature  data. 

Table  I 

Binary  Diffusion  Coef f ic idhts  as  Function  of  Temperature . 
Underlined  values  are  extrapolated;  CHgO-Og  values  cora- 
puted  from  theory;  pressure  =  1  atmosphere*? _ 

_ _ _ D(cm  /sec) _ 


T(°K) 

CH4-°2 

CM 

o 

1 

CM 

o 

o 

H20-°2 

*V°2 

C0-02 

CH2°-°2 

300 

0.226 

0. 161 

0.288 

0.821 

0.224 

0.180 

500 

0.581 

0.419 

0.692 

2.09 

0.542 

0.422 

700 

1  .  05 

0.767 

1.22 

3.76 

0.956 

0.739 

900 

1.63 

1.19 

1.85 

5.79 

1.45 

1.12 

1100 

2.29 

1.68 

2.58 

8.14 

2.03 

1.57 

1300 

3.03 

2.23 

3.41 

10.8 

2.68 

2.07 

1500 

3.86 

2.84 

4.32 

13  7 

3.40 

2.63 

1700 

4.75 

3.50 

5.32 

16.9 

4.18 

3.24 

1900 

5.72 

4.21 

6.40 

20.3 

5 .03 

3.91 

2100 

6.76 

4. 98 

7.55 

24.0 

5.94 

4.61 

The  mass  flux  fractions  ,  having  all  the  diffusion 
velocityes  V\  at  each  0.02  cm  interval  throughout  the  flame, 
were  computed  from  Eq .  (3).  These  are  shown  as  the  broken 
curves  in  Eigs  .  1  through  7.  The  great  importance  of  dif¬ 
fusion  is  obvious  in  all  cases.  It  will  be  noted  that  some 
of  the  curves  (water,  carbon  dioxide,  carbon  monoxide, 
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and  hydrogen)  have  regions  of  negative  values,  which  would 
indicate  a  net  upstream  flux  of  these  species.  While  there 
does  not  appear  to  be  any  a_  priori  reason  why  this  is  not 
possible,  the  negative  gradient  dG^/dz  which  exists  up¬ 
stream  of  these  minima  also  implies  that  these  species  are 
undergoing  a  net  disappearance  due  to  chemical  reaction. 

Since  it  is  difficult  to  see  how  this  could  happen,  it 
seems  likely  that  the  negative  regions  are  not  real  and 
should  be  regarded  as  errors  in  the  over-all  treatment.  By 

the  same  token,  the  positive  maximum  in  the  Gn  probably 

2 

is  not  real  since  it  is  a  direct  consequence  of  the  negative 
G  values  for  the  other-  species  by  the  normalization  proce¬ 
dure  used  to  obtain  the  oxygen  diffusion  velocity. 

As  in  the  previous  propane-air  flame  (Ref.  4)  it  is 
important  to  examine  the  data  in  terms  of  conservation  of 
atomic  species.  If  n.  is  the  number  of  atoms  of  a  particular 
element  in  a  molecule  of  species  i,  it  was  shown  that  the 

relation1 

I 

if  .  ' 

In,  G./M.  =  constant  (6) 

jj;  1  1  1 

must  hold  throughout  the  flame,  since  atoms  are  neither 
created  nor  destroyed.  In  the  methane-oxygen  flame,  carbon 
and  hydrogen  are  the  two  elements  whose  conservation  is  to 
be  examined,  the  oxygen  being  present  in  such  excess  in  the 
form  of  0,.,  that  its  constancy  is  practically  assured.  Re¬ 
sults  of  applying  Eq.  (6)  for  carbon  and  hydrogen  to  the 
G.  data  given  in  Figs.  1  through  6  are  shown  in  Fig.  8. 

The  data  are  given  as  per  cent  deviation  from  the  initial 
value  of  the  summation,  that  is,  the  cold  gas  composition 
upstream  of  the  screen.  (Because  of  diffusion,  the  probe 
composition  just  downstream  of  the  screen  at  z  =  0.03  cm 
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Fig.  8  CONSERVATION  OF  CARBON  AND  HYDROGEN  (ATOM  BALANCES)  THROUGH  FLAME  ZONE 
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is  not  quite  equal  to  the  true  cold  gas  composition.)  The 
balances  are  considered  quite  good, and  are  definitely  better 
than  those  in  the  earlier  propane-air  work.  The  worst  de¬ 
viation  is  -12  per  cent  for  carbon  and  -9  per  cent  for  hy¬ 
drogen,  while  in  most  of  the  flame  the  balances  are  con¬ 
siderably  better  (average  deviations  through  the  flame  are 
about  -3  per  cent).  While  the  atom  balance  is  not  a  very 
sensitive  test  oi  the  sampling  data  and  the  diffusion 
corrections,  the  fact  that  it  is  reasohably  good  certainly 
indicates  there  are  no  large  errors  in  the  data  as  a  whole. 

Conservation  of  Energy 

It  is  of  interest  next  to  examine  the  data  in  terms 
of  energy  considerations.  Since  the  flame  is  an  essentially 
constant  pressure  system,  the  convenient  energy  variable  to 
use  is  the  enthalpy,.  The  energy  continuity  equation  appro¬ 
priate  to  a  quasi-one-dimens ional  flow  with  area  change  such 
as  is  found  in  the  flat  flame  must  first  be  derived. 

The  basic  energy  equation  for  a  steady-state  system 
with  no  viscosity,  radiation,  or  external  forces  may  be 
written 

y.  (p  v  H  +  q  )  =  O  ,  (7) 

where  v  is  the  mass  average  velocity  vector,  H  is  the 
specific  enthalpy  of  the  mixture,  and  q  is  the  heat  flux 
vector  given  by  (page  717  of  Ref.  10) 

t  =  -  l\7T  +  2  N  H  t  .  (8) 

j  -L  J-  X 
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if  the  inverse  thermal  diffusion  (Dufour)  effect  is  neg¬ 
lected.  X  is  the  mixture  thermal  conductivity,  ^  is  the 
diffusion  velocity  vector,  and  H.  the  absolute  molar  en- 

1  A 

thalpy  of  species  i.  The  enthalpy  per  unit  mass  H  is  re¬ 
lated  to  the  by 


A 

H 


=  X  x,  H.  /  £  x.  M.  . 
H  ±  U  ,  11 


(9) 


Making  use  of  th^  theorem  (page  815  of  Ref.  10)  which 
states  that.,  if  F  is  a  vector  which  is  tangent  to  the  sur¬ 
face  f(x,y)  =  z,  then 

y<  V-  ^ dxdy  =  ^  Az dxdy » 

Eq .  (7)  may  be  converted  to 


+  q) 


dxdy 


constant , 


(10) 


where  v  and  q  are  the  z-components  (that  is,  normal  to  the 
flame  front)  of  their  respective  vectors.  This  may  be 
written  as 

A 

pvaH+qa=  constant,  (11) 

An 

the  quantities  p,  v,  H,  and  q  now  representing  suitably 
averaged  values  across  the  streamtube  area  a.  Combining 
the  z-component  of  Eq.<8)  with  Eq.  (11),  and  using  the 
quasi-one-dimensional.  over-all  continuity  relation  in  the 
form 


P0  vo  =  p  v  A  ’  (12) 

where  A  =  a/a  and  the  subscript  "  "  refers  to  the  cold 
o  o 

boundary  (that  is,  the  screen),  Eq.(ll)  becomes  finally 

Po  vQ  H  -  A  x(dT/dz)  +  A  S  N.  H.  V.  =  pQ  vq  $Q, 

1  (13) 
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assuming  that  gradients  in  temperature  and  composition  are 
essentially  zero  at  the  screen. 

This  equation  represents  the  enthalpy  conservation 
through  the  flame  zone  in  terms  of  the  flux  of  enthalpy 
per  unit  of  inlet  area  at  the  screen.  The  first  term  on 
the  right  side  is  the  flux  due  to  convection,  the  second 
is  that  due  to  conduction,  and  the  third  that  due  to  dif¬ 
fusion.  These  terms  were  evaluated  for  the  experimental 
methane-oxygen  data.  The  specific  enthalpy  required  was 
obtained  from  Eq .  (9)  using  tabulated  molar  enthalpies 
from  standard  sources  (Refs.  11  and  12).  The  enthalpies 
for  formaldehyde  were  obtained  from  the  heat  capacity  data 
of  Stevenson  and  Beach  (Ref.  13).  The  thermal  conductivity 
was  taken  to  be  that  of (pure  oxygen,  rather  than  attempting 
the  exceedingly  laborious  computations  for  mixtures  (for 
which  all  the  necessary  data  are  not  available  anyway) . 
Because  of  the  preponderance  of  oxygen,  this  approximation 
was  adequate  (see  Appendix  C) .  The  thermal  conductivity  of 
oxygen  has  not  been  measured  up  to  the  temperature  required 
here,  so  it  was  calculated  by  the  following  semi-empirical 
procedure :  The  translational  (or  monatomic)  thermal  con¬ 

ductivity  A.  was  obtained  from  the  NBS  tabulation  (Ref. 
mon 

11)  of  viscosity  (tj)  by  means  of  the  rigorous  theoretical 
(first  approximation)  relation 

Xmon  =  (15/4) (R/M) (r)  . 

The  improved  Eucken-type  correction  of  Hirschfelder  (Ref. 
14)  was  then  applied,  that  is 

X/Xmon  =  0,469  +  0,354  (Cp/R  "  > 

to  account  for  the  effect  of  internal  degrees  of  freedom. 
The  temperature  gradient  was  obtained  by  numerical 
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differentiation  of  the  temperature  profile,  and  the  other 
quantities  in  Eq.  (13)  were  already  available. 

The  various  enthalpy  fluxes  are  shown  in  Fig.  9. 

The  conduction  term  is,  of  course,  negative,  since  it 
always  represents  a  flux  of  heat  (positive  enthalpy)  up¬ 
stream.  The  diffusion  term  would  be  expected  to  be  posi¬ 
tive  relative  to  the  inlet,  flux  as  it  is,  since  the  react- 
ants  of  high  enthalpy  diffuse  downstream  and  the  products 
of  low  enthalpy' diffuse  upstream,  both  effects  constituting 
a  net  positive  enthalpy  flux.  (It  seems  unlikely  that  the 
effects  of  intermediates,  which  diffuse  in  both  directions 
depending  on  which  side  of  their  maximum  in  concentration 
one  considers,  would  alter  this  statement.  They  have  en¬ 
thalpies  lower  than  the  reactants  and  higher  than  the  final 
products.)  Whether  the  convection  term,  that  is,  the  spe¬ 
cific  enthalpy  profile,  constitutes  a  positive  or  negative 
flux  depends  upon  the  relative  magnitudes  of  the  conduction 
and  diffusion  terms,  as  has  been  noted  by  Lewis  and  von  Elbe 
(Ref.  15).  These  terms  depend  upon  their  respective  trans¬ 
port  coefficients.  In  the  special  case  (the  theoreticians' 
favorite)  that  all  the  binary  diffusion  coefficients  are 
equal  to  each  other  and  to  the  thermal  diffusivity  A. jv  p 
(all  Lewis  numbers  unity),  the  conduction  and  diffusion 
fluxes  cancel  and  the  convection  flux  is  constant  through 
the  flame.  This  has  been  discussed  recently  by  Hirschfelder 
(Ref..  16).  In  the  data  of  Fig.  9,  the  convection  term  shows 
a  small  positive  maximum,  but  the  present  authors  attach  no 
general  significance  to  this.  The  sum  of  the  three  fluxes 
shown  as  the  dotted  curve  would  ideally  be  constant  through¬ 
out  the  flame.  The  fact  that  the  sum  shows  a  maximum  greater 
than  the  convection  term  itself  indicates  that  the  sign  of 
the  latter  is  not  reliably  established  in  this  case. 
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APPENDIX  A 


Binary  Mixture  Approximation  for  Diffusion 


The  use  of  rigorous,  multicomponent  diffusion  coeffi¬ 
cients  in  a  mixture  even  as  relatively  simple  as  dealt  with 
here  is  completely  impractical,  a.nd  {.fortunately)  not  neces- 
ary  either.  For  any  species  i  present  in  small  concentra¬ 
tion  in  a  mixture,  it  can  be  shown  (Ref.  17)  that  ,  the 
effective  diffusion  coefficient  of  species  i  is  closely 
approximated  by 


l-mix 


1  -  x.  M./  M 
_ x  1 

£  x  ,/D  .  . 

in  3  1J 


<A1) 


This  relation  has  been  experimentally  verified  (Refs.  18  and 
19) .  If,  in  addition,  one  species  k  is  present  in  such 
excess  that  all  the  other  components  may  be  regarded  as 
traces,  then  the  above  simplifies  still  further  to 


D. 


l  -mix 


=  D. 


Dik  » 


(A2) 


which  is  the  approximation  used  in  Eq .  (4).  This  is  cer¬ 
tainly  justified.  In  the  flame  considered  here,  .in  no  case 
did  the  diffusion  coefficient  D.  „n  differ  from  the  D. 
calculated  from  Eq.  (Al)  by  more  than  3  per  cent,  which  is 
about  the  experimental  error  of  the  diffusion  coefficient 
measurements . 
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APPENDIX  B 

Neglect  of  Thermal  Diffusion 

For  any  tract  component,  the  expression  for  the 
diffusion  velocity  including  thermal  diffusion  may  be 
written  (Ref.  10) 


V. 

1 


D. 

_i 

x  . 


dx . 


k„ 


i  T  dT 

dz"  -  "T  dz 


(Bl) 


1  i- 

where  the  positive  sign  holds  if  the  trace  i  is  heavier 
than  the  main  component,  and  the  negative  sign  if  the 
trace  is  lighter, 
by 


k^,  is  the  thermal  diffusion  ratio  given 


k„  =  +  D.  /  NM.  D . 
V  —  i  li 


(B2) 


D. 

1 


being  the  trace  thermal  diffusion  coefficient,  and  N 

T 

the  total  molar  density.  In  most  cases,  D^  itself  is 
positive  if  the  trace  is  heavier  and  negative  if  it  is 
lighter  than  (in  our  case)  oxygen,  so  that  kT  is  generally 
positive.  Therefore,  the  common  situation  is  for  a  heavier 
trace  species  to  thermally  diffuse  upstream  and  a  lighter 
species  downstream. 


The  thermal  diffusion  correction  to  the  diffusion 
velocity  was  carried  out  for  the  case  of  H^O,  for  which  it 
probably  would  be  as  large  as  for  any  of  the  major  species. 
The:  necessary  thermal  diffusion  data  were  computed  from  our 
concentration  diffusion  coefficients  and  viscosities  as  out¬ 
lined  by  Amdur  and  Mason  (Ref,  20).  At  most,  the  thermal 
diffusion  contribution  to  V„  „  was  about  10  per  cent.  Thus, 
while  the  effect  is  not  always  as  small  as  one  would  like, 
the  error  caused  by  its  neglect  is  at  worst  probably  no 
larger  than  others,  and  the  simplification  this  neglect 
brings  to  an  already  complex  analysis  is  most  welcome. 
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APPENDIX  C 

Pure  Oxygen  Approximation  for  Thermal  Conductivity 


The  use  of  \q  instead  of  the  true  mixture  thermal 
conductivity  was  occasioned  primarily  (as  with  the  other 
approximations)  by  the  great  simplification  it  afforded. 

The  actual  flame  mixture  was  composed  throughout  of 
roughly  80  per  cent  oxygen,  with  the  remaining  20  per  cent 
mostly  lighter  components  (except  for  COg)  which  might  be 
expected  to  increase  the  thermal  conductivity  over  that  of 
pure  oxygen.  To  get  some  idea  of  the  error  introduced  here, 
an  approximate  calculation  (Ref.  21)  was  made  of  \  for  a 
mixture  composed  of  85  per  cent  oxygen  and  15  per  cent  water 
at  1500°K.  This  gave  a  value  about  15  per  cent  higher  than 
the  pure  oxygen  value.  Such  an  increase  in  A.  would  improve 
the  enthalpy  balance  of  Fig.  9  somewhat;. 
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WITH  A  DEFINITELY  RELATED  GOVERNMENT  PROCUREMENT  OPERATION, 
THE  U.  S.  GOVERNMENT  THEREBY  INCURS  NO  RESPONSIBILITY,  NOR  ANY 
OBLIGATION  WHATSOEVER;  AND  THE  FACT  THAT  THE  GOVERNMENT  MAY 
HAVE  FORMULATED,  FURNISHED,  OR  IN  ANY  WAY  SUPPLIED  THE  SAID 
DRAWINGS,  SPECIFICATIONS,  OR  OTHER  DATA  IS  NOT  TO  BE  REGARDED  BY 
IMPLICATION  OR  OTHERWISE  AS  IN  ANY  MANNER  LICENSING  THE  HOLDER 
OR  ANY  OTHER  PERSON  OR  CORPORATION,  OR  CONVEYING  ANY  RIGHTS  OR 
PERMISSION  TO  MANUFACTURE,  USE  OR  SELL  ANY  PATENTED  INVENTION 
THAT  MAY  IN  ANY  WAY  BE  RELATED  THERETO. 


